High speed imaging of traveling waves in a granular material during silo discharge 
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We report experimental observations of sound waves in a granular material during resonant silo 
discharge called silo music. The grain motion was tracked by high speed imaging while the resonance 
of the silo was detected by accelerometers and acoustic methods. The grains do not oscillate in phase 
at neighboring vertical locations, but information propagates upward in this system in the form of 
sound waves. We show that the wave velocity is not constant throughout the silo, but considerably 
increases towards the lower end of the system, suggesting increased pressure in this region, where 
the flow changes from cylindrical to converging flow. In the upper part of the silo the wave velocity 
matches the sound velocity measured in the same material when standing (in the absence of flow). 
Grain oscillations show a stick-slip character only in the upper part of the silo. 



When a granular material is slowly discharged from a 
vertically aligned container (tube) one often observes vi- 
brations accompanied by a booming low frequency (typ- 
ically 50 — 200 Hz) sound. This phenomenon can be ob- 
served in large silos and is often referred to as silo music 
or silo honking, but it can also easily be reproduced in 
smaller scale laboratory experiments Una- Silo quaking 
can also be observed during discharge [9|, [l0( where short 
sharp bursts of motion are separated by longer periods 
of no movement. 

Although resonances during silo discharge are long 
known in industry and have recently been investigated 
thoroughly (D-Il] , the phenomenon is not yet fully under- 
stood. While one explanation is based on stick-slip mo- 
tion of the grains near the silo walls 4 6], other authors 
argue that the phenomenon originates from a dynamic 
interaction between the silo structure and the flowing 
material, where the source of the vibration is the transi- 
tion zone between cylindrical and converging flow in the 
lower part of the system near the outlet [l|-yl • Recent lab- 
oratory experiments carefully investigated the frequency 
spectrum of the vibrations by piezoelectric accelerome- 
ters and microphones (Tj-Q and related these observations 
of the calculated eigenfrequencies of the silo structure Q 
and the time evolution of the density of the material de- 
tected by a non-invasive method (Electrical Capacitance 
Tomograph y) [Q , Q ■ Optical imaging of silo quaking was 
reported in [ljj using a frame rate of 8.3 fps. Here the 
predominant frequency of the bursts was about 1 Hz, 
which was significantly lower (by a factor of 2 or 4) than 
the frequency of the silo wall vibrations indicating no res- 
onance in this case. Moreover, grain oscillations were in 
phase over the largest part of the silo wall, so grain mo- 
tion could not be linked to wave propagation. This con- 
tradicts other measurements on silo quaking where 
the quake propagated upwards in the silo with a well de- 
fined velocity and growing amplitude, as detected by ac- 
celerometers embedded into the granular material. Very 
recent observations on resonant pipe flow (the pipe 



is continously fed from above) detected upward propa- 
gating waves using piezoelectric sensors, and claim con- 
stant wave speed in the entire tube, meaning, that the 
transition zone near the outlet has no special role in the 
formation of the resonance. This contradicts other mea- 
surements where manipulating the geometry of the 
transition zone allows for very effective reduction of the 
resonance. 

In this Brief Report we present direct observations of 
the grain motion during resonant silo discharge by means 
of high speed imaging using a frame rate of 3000 fps. 
This non-invasive technique allows us to better charac- 
terize the spatial evolution of the system. The aim of this 
study is to measure whether the grain oscillations show 
a stick-slip character near the walls, and if so, in which 
part of the silo? Are there regions in the silo, where the 
grains oscillate coherently in neighboring vertical posi- 
tions, or not? What is the speed of wave propagation in 
the granular material during resonance? 

In the present experiments the silo flow was realized in 
a vertically aligned glass tube of length 150 cm (see Fig. 
[TJa)). Two tubes were used with the inner diameters of 
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FIG. 1: (color online) (a) Schematic illustration of the exper- 
imental setup and (b)-(c) microscopic images of the materials 
used: (b) glass beads with mean diameter of d = 0.15 ± 0.05 
mm, (c) copper particles with d = 0.15 ± 0.05 mm. 
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D = 3.55 cm and 2.3 cm. We present results for two ma- 
terials (i) nearly spherical copper particles with average 
grain diameter of d = 0.15 ± 0.05 mm and (ii) spherical 
glass beads with d = 0.15 ±0.05 mm (see Figs. [Hb)-(c)). 
Qualitatively similar results were obtained with quartz 
sand, granite and corundum of similar sizes, but in those 
cases the resonance was less pronounced. For the de- 
tection of the grain motion we used high speed cameras 
(Mikrotron CAMMC1310 and 1362) with a frame rate of 
3000 fps. The resonance was also detected by traditional 
methods: piezoelectric accelerometers mounted on the 
glass tube (for horizontal vibration), on the support of 
the tube (for vertical vibration), and a condenser micro- 
phone with large diaphragm (IMG Stageline ECM-140). 
The silo discharge speed was varied by changing the di- 
ameter of the outlet D out at the bottom of the tube. As 
the main features of the phenomenon we are focusing on 
do not depend on the precise value of D out , for conve- 
nience a constant value D out = 8 mm was chosen for this 
study. 

First we present the general features of the resonance 
obtained by piezoelectric accelerometers and the micro- 
phone. In one experiment the above signals have been 
recorded during the initially full hopper discharges. This 





500 r 




450 - 




400 - 




350 - 


N 


300 - 


I 


250 - 




200 - 




150 - 




100 - 




50 - 




500 r 




450 - 




400 - 




350 - 


77 


300 - 


I 


250 - 




200 - 




150 - 




100 - 




50 - 



(a) 



1 



(b) 



0.8 






CD 


0.6 


tud 




Q_ 


0.4 


£ 

CD 


0.2 









1 




0.8 






CD 




"O 


0.6 


3 








£ 


0.4 


CD 



0.2 





10 20 30 40 50 60 70 80 90 
t[s] 



FIG. 2: (color online) Time evolution of the DFT power spec- 
trum of the signals obtained by the (a) accelerometer mounted 
on the tube, (b) microphone. Data obtained with copper 
particles in the tube with the diameter of D = 3.55 cm. The 
dashed line corresponds to the frequency of the standing wave 
(first harmonic) in the air column above the material. 

takes about 80 s and 120 s for the tube with D = 2.3 cm 
and D = 3.55 cm, respectively. The power spectrum of 
these signals has been calculated using discrete Fourier 



transform (DFT) with measurement windows of 0.5 s and 
is shown as a function of time in Figs. Eta) and (b). As 
it is seen in Fig. (J^a) at the beginning of the process the 
system vibrates with a fundamental frequency of about 
60 Hz, and this frequency slightly increases as the silo 
discharges. At t = 44 s there is a quick transition char- 
acterized by a jump in the dominant frequency, but the 
trend of a slowly increasing frequency of the dominant 
mode persists until the resonance disappears at about 
t = 78 s. This picture is fully coherent with other ob- 
servations [H-ll] and can be explained as follows. The 
granular flow excites the eigenmodes of mechanical vibra- 
tion of the tube itself and the frequency of these modes 
is slightly increasing in time as the material is flowing 
out of the system and the weight of the oscillating body 
decreases. The time evolution of the system allows for 
changing the actually excited eigenmode, for which we 
have an example at t = 44 s where we see a jump in the 
dominant resonant frequency. 

The signal picked up by the microphone shows the 
same characteristics as that of the accelerometer, but in 
addition to the eigenmode with slightly increasing fre- 
quency we find a resonance with decreasing frequency. 
This is the sign of standing sound waves in the air col- 
umn in the tube above the granular material. As the 
length I of this column is increasing with time (as the 
silo discharges) the frequency of the dominant mode is 
decreasing. The amplitude of this signal gets strong 
when it resonates with the "driving" frequency. The 
dashed line corresponds to the frequency of the stand- 
ing wave (first harmonic) according to the formula of 

/(*) = W) = im = i ■ v.t+ .3D where v * is the vclocit y 

by which the top surface of the granular bed is sinking, 
c is the velocity of sound in air and 0.31? is the length 
correction for the open end of the tube [l3| . The value 
of v s has been measured in an additional experiment, 
where the whole process was recorded optically and the 
position of the top surface was determined by digital im- 
age analysis. Several measurements confirmed, that v s is 
nearly independent of time except near the very end of 
the process where v s slightly increased. 

In the following we will focus on the data obtained 
with high speed imaging. First we show how the grains 
move at a selected height x. For this we took images of a 
vertical stripe of 1 x 0.1 cm and determined the displace- 
ment of the particles by correlating subsequent images. 
This procedure was repeated at various heights. In Fig. 
|3] three examples are shown, where Fig. shows the 
position s of the particles as a function of time, while (b) 
and (c) show the corresponding velocity v and accelera- 
tion a. The three measurements correspond to x = 12 
cm, x = 67 cm and x = 127 cm measured from the bot- 
tom of the tube. As it is seen, the grains move downwards 
with oscillating velocity. At the top and middle of the 
tube during each period they come to a full stop, remain 
standing for a considerable portion of the period and then 
accelerate downwards. Their acceleration almost reaches 
the value corresponding to free fall (g = —9.81 m/s 2 ) 
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near the top of the tube, but is lower in the middle of 
the tube. At the bottom of the tube one can still see 
oscillations in the grain velocity but the grains do not 




0.04 0.06 
t[s] 

FIG. 3: (color online) (a) Grain position, (b) velocity and 
(c) acceleration as a function of time taken at three heights 
x — 12 cm (x), x = 67 (A) cm and x = 127 (o) cm measured 
from the bottom of the tube. Data taken with copper particles 
in the tube with the diameter of D = 3.55 cm. 

come to a full stop and their acceleration does not come 
close to g. The average grain velocity near the walls v av 
(i.e. the slope of the curves in Fig. 02a)) is smaller than 
the velocity v s by which the top surface of the granu- 
lar bed is sinking. The corresponding three values are 
Vav/v s = 0.54, 0.59 and 0.89 for x — 12 cm, 67 cm and 
127 cm, respectively. This means that no perfect plug 
flow develops, i.e. the grains move slower near the walls 
than in the center of the tube even at the uppermost part 
of the material. 

In order to measure whether or not the grains at dif- 
ferent heights oscillate in phase with each other, a longer 
vertical stripe of 13.2 x 0.5 cm was monitored. Two 2.2 
cm long segments were analyzed at the upper and lower 
parts of this image. As it is seen in Figs. Hfa)-(b) there is 
a phase shift between these signals (indicated by vertical 
lines) implying that the grains in neighboring segments 
do not oscillate in phase with each other, but waves are 
propagating upwards in the system. The phase shift be- 
tween the neighboring peaks in Fig. H{a) corresponds to 
a wave velocity of U = 54 m/s, while in the lower part of 
the tube (Fig. QJb)) a larger wave velocity (smaller phase 
shift) is observed. The velocity U of the upward prop- 
agating waves has been measured at 6 different vertical 
positions. The resulting U (x) curves are shown in figs. 
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FIG. 4: (color online) Grain velocity v as a function of time 
at two nearby locations at (a) x = 102 cm (*) and x — 113 
cm (o) and (b) x = 22 cm (*) and x = 33 cm (o). The 
time lag between the curves is a measure of the corresponding 
sound velocity, (c)-(d) Velocity of sound waves as a function 
of location x in the tube with (c) D — 2.3 cm and (d) D = 3.55 
cm. The curves correspond to different stages of the process, 
when the top of the material was at x = 145 cm (o), x = 110 
cm (V), x = 11 cm (o) and x = 47 cm (A). Data for copper 
with open symbols, data for glass beads with solid symbols. 



|4jc) and (d) for the two tubes with D = 2.3 cm and 3.55 
cm, respectively. This measurement has been repeated 
at four different stages of the process, where the top of 
the material was at x = 145 cm, x = 110 cm, x = 77 
cm and x = 47 cm. The wave velocity U appears to be 
constant in the upper part of the tube. In the smaller 
tube (D = 2.3 cm) we get U — 44 ±6 m/s for glass beads 
and U = 54 ± 6 m/s for the copper particles, while in 
the larger tube (D = 3.55 cm) U = 65 ± 6 m/s for both 
materials. These values are similar to those previously 
observed for the velocity of upward propagating pulses 
detected during "silo quaking" There single rarefac- 
tion waves (quakes) propagated upwards in the system 
with a velocity of 74 ± 15 m/s and 135 ± 20 m/s in the 
two Pcrspex tubes with the diameters of 9 cm and 30 cm 
respectively. The fact that we find larger value of U in 
the larger tube also agrees with the above data obtained 
for quakes and is in accordance with very recent observa- 
tions on continous pipe flow, where the velocity of sound 
waves was demonstrated to increase with increasing tube 
diameter fl2j . Our more detailed spatial characterization 
however revealed (Figs. IUc)-(d)), that the wave veloc- 
ity is not constant throughout the tube but it increases 
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towards the bottom of the tube. The larger value of U 
measured at the lower part of the tube means a smaller 
phase lag between the oscillation of neighboring particles, 
which, in accordance with other results [l|, Q , shows that 
the resonance originates from that part of the tube where 
the flow changes from cylindrical to converging flow (at 
x rj 4D) leading to strong density oscillations. 

An additional experiment was done to measure the ve- 
locity of sound waves in the standing material (without 
flow) using the smaller tube (D = 2.3 cm). In this mea- 
surement we placed 2 piezoelectric accelerometers in the 
granular material at a distance dx above each other, and 
perturbed the granular bed at the lower end of the tube 
with a pneumatically generated single pulse. This was 
done using a smaller (flexible) tube with a rubber mem- 
brane at the end, which was placed in the lower end of 
the silo. We measured the time dt needed for the pulse 
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FIG. 5: (color online) The elapsed time dt during a sin- 
gle pulse traveled between two piezoelectric accelerometers 
placed in the static material at a distance dx for the case of 
copper (A) and glass beads (o) in the tube with D = 2.33 cm. 
The slope of the linear fits (shown with dashed lines) gives U a 
the velocity of sound waves in the standing material. 

to travel between the two accelerometers as a function of 
dx. The slope of the resulting curves (presented in Fig. 
[5]) gives for the speed of sound waves in the standing ma- 
terial U s — 53 ± 5 m/s for glass beads and c7 s = 48 ± 5 
m/s for the copper particles. These values agree nicely 
with the ones discussed above observed in the same tube 
during flow (resonance) and are similar to other exper- 
imentally obtained values for the sound speed in glass 
beads (60 - 240 m/s 0; 40 - 60 m/s 0) or sand (40 
m/s |l5|). or even using much softer photoelastic parti- 
cles (100 m/s [13 )■ These sound speed measurements 



were done in the absence of strong external load and 
also revealed, that the result can significantly depend on 
the measurement method More importantly, recent 
investigations uncovered, that in samples under exter- 
nal uniaxial load, the velocity of compressional waves is 
sensitive to stress induced anisotropy, leading to larger 
U perpendicular to the applied load than parallel to it 
[19|. Also, U increases more gradually with pressure P 
(U oc P 1 / 4 ) at lower pressures, compared to the generally 
observed U oc P 1 / 6 dependence at higher P [T^]. Both of 
these observations [H[[li| indicate, that the larger wave 
velocity measured during resonant flow in the lower part 
of the tube suggests increased radial pressure in this re- 
gion, where the velocity gradient changes due to the prox- 
imity of the outlet. 

In summary, we studied a resonance phenomenon 
called "silo music" during silo discharge. We have vi- 
sualized the oscillations in grain motion with high speed 
imaging, and detected that grain oscillations are not co- 
herent when compared at different heights, but there is a 
phase lag. Information propagated upward in the system 
in the form of sound waves. One of our major findings 
is, that the wave velocity is not constant throughout the 
silo, but considerably increases towards the lower end of 
the system (corresponding to a more coherent motion of 
the grains). This underlines the importance of the tran- 
sition zone in the lower end of the silo, where the flow 
changes from cylindrical to converging flow and is in ac- 
cordance with those observations [JQ which suggest that 
the source of the resonance phenomenon is the generation 
of strong stress oscillations in this region. As shown, the 
amplitude of the grain oscillations increases with height 
and leads to stick-slip motion in the upper part of the 
silo. This supports the argument that appropriate fric- 
tion with the walls is crucial for the amplification of the 
waves 0, [TH, HU HH and could be a necessity for strong 
resonance as argued by other authors d-Q , which is also 
supported by the fact that increasing the roughness of 
the silo walls (i.e. suppressing stick-slip) effectively re- 
duces the resonance @, Q . Finally, we have also shown 
that the velocity of the waves in the upper part of the 
silo matches the sound velocity measured in the same 
material when standing (in the absence of flow) . 

The authors are thankful for discussions with Kesava 
Rao, who motivated us to study this problem. We thank 
A. Konya for critical reading of this manuscript. 
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